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Abstract 

In relativistic heavy-ion collisions, the strong Lorentz-contracted electromag- 
netic fields are capable of producing copious numbers of lepton pairs through 
the two-photon mechanism. Monte Carlo techniques have been developed 
that allow the exact calculation of production by this mechanism when a 
semi-classical approximation is made to the motion of the two ions. Here, 
we develop a hybrid Monte Carlo technique that enables us to calculate the 
impact parameter dependence of the two-photon mechanism for lepton-pair 
production, and by using this result, we obtain the probability distribution 
for multiple-pair production as a function of impact parameter. Computa- 
tions are performed for S-|-Au and Pb-|-Pb systems at 200 A GeV and 160 A 
GeV, respectively. We also compare our results with the equivalent photon 
approximation and elucidate the differences. 
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I. INTRODUCTION 



The proposed new colliding-beam accelerators - designed to investigate nuclear matter 
at high temperatures and densities by accelerating highly charged heavy ions at fixed-target 
energies per nucleon up to 20 TeV - have motivated great interest during the last decade 
concerning possible new electromagnetic phenomena. When heavy ions collide at relativistic 
velocities, the Lorentz-contracted electromagnetic fields in the space-time region near the 
collision are sufficiently intense to produce large numbers of electron-positron pairs, muon 
pairs, vector bosons, and possibly the yet-unconfirmed Higgs boson. All these processes 
occur at nearly atomic distance scales The phenomena involved are pervasive, im- 

pinging upon atomic, nuclear, and particle physics. The electromagnetic fields associated 
with these collisions are intense as they are proportional to the Lorentz factor 7, which 
is approximately the beam kinetic-energy in units of GeV per nucleon, and the charge of 
the ion Z. These parameters, together with the impact parameter, 6, determine the fields 
available in a collision. The feature of electromagnetic particle production by heavy ions, 
which piqued early studies, is the Z'^ enhancement through the coherent action of all the 
constituent charges of the colliding partners. The pulsed fields are strongly time dependent 
with a width of 6/7 and, thus, contain large Fourier components which give rise to rela- 
tively large particle production probabilities. Coherent particle production is most clearly 
distinguished in peripheral collisions. 

One of the most interesting of these processes, from the perspective of both fundamental 
and practical importance, is the sparking of the vacuum to produce electron-positron pairs. 
Much study has been devoted to the problem of electron-positron pair production during 
recent years in anticipation of new experimental opportunities at the Relativistic Heavy-Ion 
Collider (RHIC), currently under construction at Brookhaven National Laboratory, as the 
immense predicted fiuxes of electrons cannot be ignored in detector or accelerator design. 
This facility will provide extremely relativistic colliding beams of fully stripped ions as heavy 
as gold, fully exposing the large charge of the atomic nucleus. 

The primary physics goal of the RHIC project is the creation and study of a so-called 
quark-gluon plasma. This unique form of matter is expected to be formed in the central, or 
near-central, collisions of heavy ions at extreme relativistic energies [^. The thermodynamic 
conditions attained in these central collisions are expected to be such that the constituent 
quarks and gluons of baryons and mesons become deconfined in a new and short-lived plasma 
state. Electron- and muon-pair production from hadronic interactions have been widely 
discussed as a possible tool to help probe the formation and the decay of the quark-gluon 
plasma phase of matter [|,^ . In the conditions of such central collisions, lepton-hadron final- 
state interactions are usually small, and, hence, the leptons carry direct information on the 
space-time region of creation. However, suggestions by several authors indicate that other 
sources of lepton pairs might possibly mask the leptonic signals originating from the plasma 
phase • Electromagnetic production from the vacuum of single- and multiple-lepton 
pairs is a major contribution to this physical background pil|JI[] and, therefore, must be 
understood in detail. In addition, two very abundant electromagnetic processes constitute 
the primary limitation to the lifetime of stored beams at RHIC. One is a nuclear decay 
following the electromagnetic excitation of the giant dipole resonance, and the second is 
the creation of an electron-positron pair accompanied by the capture of the electron in an 
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atomic bound state of a participant heavy ion. Both processes result in a change in the 
charge-to-mass ratio of the ion in the storage ring causing it to be deflected out of the beam. 

From a fundamental perspective, relativistic heavy-ion collisions provide an opportunity 
to study nonperturbative quantum electrodynamics (QED) in an entirely new and contin- 
uously varying energy regime, using an interaction which is completely known due to the 
combination of very high collision energies and electric charges. The coherent, electromag- 
netic production of electron-positron pairs using heavy ions is fundamentally different from 
other production mechanisms using light particles at high energies, since in the former, the 
coupling constant is strongly enhanced due to the large charge. 

Nonetheless, perturbative methods are useful in studying the problem of lepton-pair 
production, and reliable lowest-order perturbative calculations have been used as input into 
design models for RHIC [p|,p!l|. Beginning with quantum-field theory, we derive a classical- 
field method based on Feynman perturbation theory in the limit that the momentum transfer 
of the photons is much smaller than the momentum carried by the nuclei • This is 
manifestly true for the coherent production of particles through two-photon processes in 
relativistic heavy- ion collisions. These lowest-order diagrams (see Fig. 0), coupling lepton 
fields to classical electromagnetic fields, have been evaluated exactly using Monte Carlo HJTSU 
and analytical techniques , and predict cross sections consistent with experiments for 

free electron-positron production performed by two independent groups using collisions of S 
+ Au and S + Pt [T^, respectively, at fixed-target energies of 200 A GeV. Furthermore, 



no indication of nonperturbative effects is observed in these experiments using the relatively 
light sulfur projectile. However, applying lowest-order perturbation theory to the production 
of electron-positron pairs with heavy ions at high energies and small impact parameters 
results in probabilities and cross sections which violate various theoretical bounds such as 
unitarity pT|J^. It is therefore clear that low-order perturbative calculations alone are not 
adequate for smaller impact parameters at RHIC energies, as higher-order damping effects 
must be included (e.g. the creation of multiple pairs [^). Heavy ions have been used to 
produce electron-positron pairs in collisions of U -f- Au at fixed-target energies of 0.96 A 
GeV, and the observed cross sections are not in agreement with the low-order calculations 



performed at these energies [16]. New experiments using a variety of targets are scheduled 



for November, 1994 [17| using lead beams at fixed-target energies of 160 A GeV, and for 



1996 |T^ using gold beams at fixed-target energies of 12 A GeV during which electron-pair 
multiplicities and other nonperturbative features of these collisions will be observed. 

The Monte Carlo method described in Ref. [0] for evaluating the two-photon Feynman 
diagrams employs an analytic integration over the impact parameter and, thus, does not pro- 
vide information on the impact-parameter dependence of the cross section. In this present 
work, we generalize these techniques in order to calculate the impact-parameter dependence 
of the two-photon diagrams which has implications for the description of the electromagnetic 
background for hadronic interactions, the study of strong-field effects (e.g. pair multiplici- 
ties), and the study of the validity of the equivalent-photon approximation, as is discussed 
below. Various possibilities for experimental study of the impact-parameter dependence for 



free electron- and/or muon-pair production are discussed in Ref. ||23 



Concerning certain classes of detectors designed for nuclear physics experiments at RHIC, 
the total integrated cross section may not be the most relevant quantity for studies of the 
electromagnetic background. If the detector triggers on a limited range of small impact 
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parameters, the pair production probability for these impact parameters would be more 
relevant [0]. Due to the nature of the electromagnetic interaction, the particle produc- 
tion process occurs over a large range of impact parameters. The larger impact parameters 
produce weaker, more perturbative, fields, but these impact parameters are favored by the 
geometrical weight factor in the cross-section integration. An integration over all impact 
parameters also includes contributions from impact parameters where the two nuclei would 
collide. The assumption that the nuclei continue on straight-line paths for these nonperiph- 
eral impact parameters is not accurate. Thus, one would like to isolate the relatively small 
impact parameters which are still larger than the grazing impact parameter for the study of 
strong-field effects. 

Recent progress has been made in understanding the nonperturbative nature of the pro- 
duction of multiple free electron-positron pairs p^- p7|] . Using the usual assumptions of 
the classical, strong nature of the electromagnetic field generated by the heavy ions and 
omitting the final state interactions among the produced leptons, one can express the prob- 
ability for producing A^— pairs as a Poisson distribution whose mean value is the probability 
for producing a single pair in lowest-order perturbation theory. Therefore, the central in- 
gredient needed for describing the multiplicity distribution of electron-positron pairs within 
this approximation is the impact-parameter dependence of the two-photon mechanism. 

Historically, the two-photon process has been modeled through the equivalent-photon ap- 
proximation, which takes advantage of the close relation between the interaction produced by 
a relativistic charged particle and those due to incident electromagnetic waves [p8l -|36[| . In the 
equivalent-photon approximation, the equivalent-photon flux associated with a relativistic 
charged particle is obtained via a Fourier decomposition of the electromagnetic interaction 



| 37| , |38|| . Cross sections are obtained by folding the elementary, real two-photon cross section 
for the pair-production process with the equivalent-photon flux produced by each ion. The 
pair-production cross section is easier to compute using the equivalent-photon approxima- 
tion than by calculating the two-photon diagram, and the results for the total cross section 
are reasonably accurate, provided the incident-particle Lorentz factor 7 is much greater than 
one, and that the energy transferred via the photon is much less than 7. However, details of 
the differential cross sections, spectra, and impact-parameter dependence differ, especially 
when complicated numerical cuts in the coordinates are applied in order to compare with 
experiments. Among the shortcomings of this approximation is an undetermined parameter 
which corresponds to the minimum impact parameter or the maximum momentum trans- 
fer, which makes it difficult to get specific results. As such, the method loses applicability 
at impact parameters less than the Compton wavelength of the lepton pO[, which is the 
region of greatest interest for the study of nonperturbative effects. Reference PD[ 



reviews 



the application of the equivalent-photon method in relativistic heavy-ion collisions through 
1988. 

Section |1| introduces the basic two-photon approach for computing the single-pair pro- 
duction and the generalization of these results for calculating multi-pair production cross 
sections. In this section we also discuss the impact-parameter dependence of these cross- 
sections. In Section |T| we discuss the details of the numerical techniques used for obtaining 
impact-parameter-dependent cross sections. Section ^ outlines the results for S+Au and 
Pb-|-Pb heavy-ion reactions. The paper concludes with the discussion of the results in Sec- 
tion |V]. 
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II. FORMALISM 



A. Single Pair Production 

The cross section for lepton-pair production in relativistic heavy-ion colhsions via the 
two-photon process has been derived in using a semiclassical least-action principle. The 
result is equivalent to utilizing the leading-order Feynman diagram and taking the classical 
limit 0,^,0 on the relative motion of the two ions. In this formalism, the source currents 
appear as arising from Lorentz-boosted charge distributions. We follow the formalism of , 
only here we develop a technique for calculating cross sections as a function of the impact 
parameter, rather than the integral over all impact parameters, as was done previously. The 
lowest-order two-photon process is pictured in Fig. |l[ 

Nucleus 1 is the nucleus which (in the center of velocity frame) moves with a velocity 
and nucleus 2 (in Fig. |^) moves with velocity +(3, both parallel to the z-axis. Their 
trajectories are taken to be straight lines separated by an impact parameter b. Throughout 
this paper, we use a system of units with h = c = m = e = l where m is the electron's rest 
mass. The semiclassical coupling of electrons to the electromagnetic field is given by the 
Lagrangian density 

Cintix) = -^{x)j^<i!{x)A^{x) , (2.1) 
which only depends on the field variables via the classical four-potential A'^, where 

A^" = A^il) + A''{2) , (2.2) 
and in the momentum space, the nonzero components of the potential from nucleus 1 are 



^ 2 



^ (1) = -Stt ^7 exp 

A'{1) = (3A\l). (2.3) 

The potentials from nucleus 2 can be obtained from (3) by the substitutions b — b, /5 
—(3. If we assume that the heavy-ion motion can be localized along definite impact param- 
eters, we can write the total inclusive singles cr^ and pair cTp cross sections: 

a, = 1 d%Ms{h) , 

cTp = j (fhNpih) , (2.4) 

where A/'s(b) is the singles multiplicity and Np{h) is the pair multiplicity of produced elec- 
trons. In the strong field limit the multiplicity of electrons represents the mean number of 
electrons produced out of the vacuum. It has been shown that, in the perturbative hmit, the 
pair multiplicity is the same as the singles multiplicity, and the inclusive pair cross section 



IS CTp ~ 0", 



a 



jd'hY.T.\<xt\S\x^i^>\\ (2.5) 



k>Qq<0 
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where the summation over the states k is restricted to those above the Dirac sea, and the 
summation over the states q is restricted to those occupied in the Dirac sea. The transition 
matrix element for direct Feynman diagrams in Eq. ( p.5| ) has been derived in as 



2/3 J (27r) 



- exp < t 



bUW(A:,g:pj 



(2.6) 



Including both the direct and crossed Feynman diagrams, the results for the cross section, 
as a function of impact parameter, can be obtained as 



da 1 



E 



&Jo(&|P±-P±l) 



{k,q:p±)+ A^-\k, g : + q± - p±; 
(k, g : p1) + (k, g : + - p'^; 



(2.7) 



where k (q) is the momentum of the produced lepton (anti-lepton), and A^'^^k, q;p±) are 
given by 



^(+)(/c, g : px) = F(kx - p± : uJi)F{p^ - q^ : uj2)Tkg{pi_ : +(3) 
A'^'\k, q : p_l) = -F(k_L - p_L : uj2)F{p± - q± : LJi)Tkg{p± : -/5) 



(2.8) 



The quantity -F(q, uj) is the scalar part of the electromagnetic field of the moving heavy ions 
in momentum space 



F(q : u) 



up- + /5272|qp 



GE{i)Jz{q'' 



(2.9) 



where GE{q^) is the form factor of the nucleon and fz{q^) is the form factor of the nucleus. 
The frequencies uji and uj^ of the virtual photons are fixed by energy conservation at the 
vertex where the photon is absorbed 



ijJ2 



2 

E(-) - - /3(g, - K 



(2.10) 



The quantity T contains the propagator of the intermediate lepton and the matrix elements 
for the coupling of the photon to the leptons 



T,,(px:/3)=EE 



X <<)|(1 -/?«.)!<) ><«(;)|(l + /5a,)|n(;) > . (2.11) 
Here, \u^^^ > is the usual Dirac spinor and is the z-component of the Dirac matrices. 
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B. Multiple Pair Production 



In this section we will provide the basic formulas used in our calculation of multi-pair 
cross sections. Recently, a number of differing approaches have been developed to compute 
the multi-pair cross sections 0-|2^, all resulting in a Poisson distribution for the multi- 
pair impact-parameter-dependent probability function. The prevailing assumption is that 
the heavy ions do not suffer any recoil or energy loss while electron-positron pairs are pro- 
duced from the vacuum, which seems to be a reasonable assumption when the pair energy 
is compared with that of the heavy ions. One implication of the above assumption is the 
approximation of straight-line, heavy-ion trajectories for all impact parameters. Although 
this may not be correct for central impact parameters, its reliability for cross-section pre- 
dictions depends on the physical region for producing pairs. Provided that most of the 
pair production cross section stems from peripheral heavy-ion trajectories, the use of this 
assumption seems to be safe. 

It has also been recognized that in most lowest-order perturbative calculations, pair pro- 
duction probability violates unitarity for small impact parameters and extreme energies. It 
has been suggested that the summation of the classes of diagrams resulting from the inde- 
pendent pair approximation can be used to restore unitarity to the lowest-order perturbation 
theory. These have been discussed in Ref. using the sudden and quasiboson approxi- 
mations, in Ref. who obtained a similar result by a straightforward summation of the 



diagrammatically defined series in the perturbational expansion, and in Ref. [^] using a 
nonperturbative treatment and neglecting the interference terms. The resulting multi-pair 
production probability is described by a Poisson form 



where denotes A^-pair production and 



nb) = j:j:\<xrwr>\"^ (2-13) 

A;>Og<0 

is the lowest-order perturbation result for the pair-production probability. P(b) can also 
be interpreted as the average number of electrons produced out of the vacuum. The A"- 
pair cross section a^pair is obtained by integrating the A^-pair probability over the impact 
parameter h. 



<yNva.r = J d'bP^ib) , A^ = 1, 2, . . . (2.14) 
In the two-photon formalism, as discussed above, the total one-pair cross section is given by 

(T = J (fhV(h) . (2.15) 
Therefore, from the above equation, we can write V (b) 
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which can be used in Eq. ( p.l2|) to calculate A^-pair cross sections. Consequently, our task 



in this paper is to determine a well-behaved impact-parameter-dependent cross section. In 
terms of the multi-pair cross section, the total cross section for producing any number of 
pairs is given by: 

N=l 



III. NUMERICAL TECHNIQUES 

For the calculation of the total cross section, the differential cross section given by Eq. 
can be integrated over all impact parameters analytically (it produces a simple delta 
function) and yields the result for the total cross section as given in Ref. However, the 
impact-parameter-dependent cross section is numerically much more difficult and evaded 
any computations until the present. This is due to the fact that the function Jo(& |p_L — pYl) 
is a rapidly oscillating function, especially for large h. One can add an additional d?p'^ 
integration to the Monte Carlo technique given in Ref. and attempt to do the integration 
in Eq. ( |2.7| ) directly. We tried up to 100 million Monte Carlo points, and, although accuracy 
for small impact parameters was acceptable, we have not been able to obtain convergence 
for large impact parameters. 

We propose the following technique for circumventing this difficulty. We divide the 
integration according to: 



— = / dqqbJo{qb)J^{q) 
do Jo 

where J-'{q) is given by a nine-dimensional integral 

72;^ j2t^j2 



(3.1) 



TT 



JO 



2- r dkJq.d'k^d'Kd'Q 
g 



X |F 
+ F 
X |f 
+ F 



-(Q-q);u;i 

^(Q - q);^i 
^(Q + q);^i 
:(Q + q);^^i 



(27r)io 
F [-K; 002] 

F[-K;uj2]%, [k^-K;-/3]| 



kx--(Q-q);+/5 



F [-q - K; ^2] % 



kg 



-(Q + q);+/5 



F [-q - K; ^2] [k^ + q - K; 



(3.2) 



where we have changed variables to 



P± = k± - -(Q - q) 

q± = k_L 
p1 = 



Q - q) + K 



2(Q + q)- 



(3.3) 
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For a fixed value of g, tlie Monte Carlo technique of Ref. can be generalized to calculate 
the nine-dimensional integral of Eq. ( p.2|) . The details are provided in the Appendix. 

This procedure results in a function, J^{q), which is a relatively smooth function of q. 
The number of Monte Carlo points used is chosen such that there is approximately a five 
percent error in each point we generate for J^{q)- We fit a smooth function to the T{q) 
generated by the Monte Carlo integration and then perform the q integral analytically. We 
check the procedure by integrating over the impact parameter and find that we are always 
within the five percent error we set for ourselves in comparison with the total cross section 
calculated using the methods of Ref. 0] . 



IV. RESULTS 

To demonstrate the ability of the technique presented here to produce the impact pa- 
rameter dependence of the two-photon electron-positron pair production cross sections, we 
present results for two typical cases. These are for electron-positron pairs produced in S 
-|- Au collision at 200 A GeV per nucleon, and in Pb + Pb collision at 160 A GeV per 
nucleon in the fixed-target frame. These two systems are studied because they correspond 
to experiments either completed or planned for the immediate future at CERN ||T^ , [T7|| . In 
all cases we set the form factor of the nucleon equal to 1 and use a Woods-Saxon form factor 
for the nucleus. The use of regular form factors for the nucleon, or uniform or Gaussian 
form factors for the nucleus does not alter the results at these energies. 

In Fig. ^ we present the function J^{q) from both the Monte Carlo calculation and from 
our smooth fit. The smooth functions fit to J^{q) are given by, for Pb + Pb and S + Au 
cases, respectively 

^(g) = 3835.7 e-i-34°29,^ 

^(g) = 163.5 . (4.1) 

The differential pair production da/db is then calculated by integrating the smooth fit to 
J-'{q) according to Eq. ( p.l|) . The results are presented in Fig. | for S + Au and Pb -|- Pb 
cases, respectively. We see a smooth exponential fall for large b as expected. We may further 
integrate these results over impact parameter. In Fig. § we show the integrated cross section 
when integrated to a maximum impact parameter, bmax, given by 

i-OO rbmax 

(T = dqq dbb jQ{qb)T{q) 

Jo Jq 

dqbraax Jl{qbmax)^{q) (4.2) 



which as b goes to infinity reproduces the total cross section obtained in Ref. using the 
exact integration over b. Analytically one could also show that jF(g = 0) converges for large 
bmax to the total cross section of 0. This is done by taking the b^ax — ^ oo limit of Eq. 
as follows 



/ dxJi{x)T{- ) , X = qb„ 

OfYiax 

lim J^(0) / dxJi{x) 

bmax—*00 JO 



^(0) , (4.3) 



9 



where we have taken T out of the integral, since, in the range of x for which the integral 
converges, T is essentially a constant with argument very close to zero. The same result 
can also be obtained by actually setting g = in the definition of T{(i) given by Eq. (p.2D, 
and realizing that the resulting expression is the same as the total cross section expression 
given in Ref . . Numerically, we have obtained this result in Eq. ( [4.1| ) with a less than two 
percent error. 

We compare our two-photon results in Figs. ^ and ^ with results obtained using the equiv- 



alent photon approximation as expressed in |39,25]. Using this equivalent photon method. 



a simple expression for is given by the lowest-order perturbative result as 



X 



c 



h 



2 



In 



c 



2b 



HZ) , (4.4) 



where 7 is approximately the energy of one ion measured in the frame of the other ion 
(7 = 2'^1^i — 1), (5 is a constant that has a value of 0.681. A(Z) is a Coulomb modification 
which is ignored here. In this equation, for impact parameters smaller than A^, V{h) goes 
to infinity very rapidly; i.e., the equivalent-photon formula diverges for very small impact 
parameters [Q. In addition, impact parameters larger than ■jSXc/2 exceed the region of 



validity for this formula. However, for valid impact parameters, i.e. 'y6Xc/2 > b > Xc, the 
simple equivalent-photon method gives quite reasonable results although it overestimates 
the contribution to pair production from the two-photon process by approximately five 
percent in both cases considered, with the largest disagreement coming at small impact 
parameters. In Fig. ^ the equivalent-photon probabilities have been integrated in their 
region of validity, whereas the integration limits for the Monte Carlo evaluation of the 
two-photon diagrams are from zero to infinity. We note that other formulations of the 
equivalent-photon approximation may implement slightly different approximations from the 



ones used to obtain Eq. ( [4.4|) pO| , and, as a result, may produce results somewhat different 
than the simple formula used here. 

Taking advantage of our demonstrated ability to compute the impact-parameter de- 
pendence of the two-photon mechanism and the simple Poissonian form for the impact- 
parameter dependence of the multiple-pair cross section (Fig. ( |2.12| )), we compute the prob- 



ability distributions for A^-pair production as a function of impact parameter b for S + 
Au (200 A GeV) and Pb + Pb (160 A GeV) collisions. Our results show that the A^-pair 
production probability has a finite value at impact parameter 6 = and is continuous every- 
where. In Fig. ^ and Fig. |^ we see that one-pair production completely dominates both the 
S -|- Au and Pb + Pb collisions. The two- and three-pair probability distributions are small 
and, as expected, decrease rapidly for impact parameter values larger than the Compton 
wavelength, although for the Pb -|- Pb case the two-pair cross section has an appreciable 
95b value. The other values for two- and three-pair cross sections are shown in Figs. |] and 

For comparison, we also plot in Fig. ^ and Fig. ^ the multiple-pair, impact-parameter- 
dependent probability distributions calculated using the equivalent-photon approximation 
(Eq. (|4.4| )). However, the inability of the equivalent-photon approximation to accurately 
approximate the two-photon mechanism at impact parameters below one Compton wave- 
length, along with the general overestimation of the total pair-production cross section by Eq. 
([4.4|), results in the equivalent-photon approximation overestimating the N-pair multiplici- 



10 



ties. For the lead collision considered, the equivalent-photon approximation overestimates 
the 2-pair cross section by 40% and the 3-pair cross section by a factor of 2. Once again, 
the equivalent-photon approximation has been integrated here only over the range of valid 
impact parameters. 

In summary, we have found the total one-pair cross section for S + Au to be 160 barns and 
Pb + Pb to be 3569 barns by integrating the two-photon diagram and using the Poissonian 
distribution to extract the cross section for one-pair to be produced. The cross section for 
any number of pairs can be found by summing the contributions from all multi-pair cross 
sections as given by Eq. |2.17| . Another interesting quantity is the total cross section for 
multi-pairs, i.e. excluding the single-pair cross section. For the lead collision considered, 
this cross section has an appreciable value of approximately 100 barns which is about 3% of 
the total cross section. 



V. CONCLUSIONS 

We have found a technique which is capable of calculating the semiclassical two-photon 
mechanism for lepton-pair production as a function of the impact parameter. This required 
a generalization of the previously developed Monte Carlo technique in order to be able to 
handle the oscillating Bessel function which occurs in Eq. ( p. 71 ). We find that utilizing a 
Monte Carlo technique to generate the remainder of the integral, J-'{q), as a function of q, 
fitting a smooth function to this result, and then performing the final integral utilizing the 
smooth fit function, produces a workable hybrid Monte Carlo technique. 

Using the simple results obtained recently in understanding the nonperturbative nature 
of the production of multiple, free electron-positron pairs [p^-p7|, we use the calculated 



impact-parameter dependence of the two-photon diagrams as the central ingredient needed 
to calculate the multi-pair production cross sections. Our calculations show that, for the 
systems and energies considered, single-pair production is the dominant part of the total 
production cross section as expected, since, at these energies, the two-photon mechanism 
does not violate unitarity, indicating that nonperturbative effects remain relatively small. 
For the Pb + Pb case we find a substantial two-pair cross section of 95 barns, which may 
yield itself to measurement. However, the impact-parameter dependence of the two- and 
three-pair cross sections are limited to a few Compton wavelengths (6 < lOAc). For S + 
Au and Pb + Pb collisions, we have found that a large portion of the cross section of e~e^ 
pair production resides in the region b > 2R (R ~ 6.8 fm for Au, Compton wavelength 
for electron Ac=386 fm); thus, for low invariant masses we can eliminate much of the 
concern over the hadronic debris. The simple equivalent-photon approximation does a fairly 
good job of predicting the total pair-production cross section, but its utility decreases when 
multiple-pair cross sections are desired due to its inability to accurately describe the small 
impact-parameter collisions. 

It appears that for the lepton-pair production problem, classical field treatment and 
quantum field theory have complementary roles. In the classical field treatment, we assumed 
that the momentum transfer of the photons is much smaller than the momentum possesed 
by the nuclei; therefore, this is in good agreement with the quantum field theory where the 
energy and momentum conservation is exactly satisfied at all vertices of Feynman diagrams. 
In relativistic heavy-ion collisions, coherent particle production via two-photon processes 
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does allow us to make the same approximation because the nuclear form factor constrains 
the virtual photon momenta where such an approximation is particularly good. Quantum 
field theory provides guidelines for the classical field treatment, and tells us what is the 
classical counterpart, especially for nuclear systems with spin, and when such a treatment 
is vahd. Classical field treatment also produces a clear physical picture for the classical 
trajectories of heavy ions which can be tagged experimentally, and enables us to properly 
simulate the dependence of particle production in impact parameter space. 

With the advent of higher-energy heavy-ion colliders, the study of the physics of the 
two-photon process emerges as an exciting new field. The process may be used as a means 
of production of exotic particles - perhaps the study of non-perturbative effects in QED. 
In all cases, an understanding of the impact-parameter dependence of the production is 
necessary, and it is useful to be able to work in regions that are beyond the validity of the 
equivalent-photon approximation. 
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APPENDIX A: MONTE CARLO INTEGRATION 



The J^(q) function in Eq. ( |3.2| ) is of the form 

J='{q) = Fo j f{xi,X2, xg, q)dxidx2 ■ ■ ■ dxg , (Al) 
where we denote these coordinates by the eight-dimensional vector 

^ = {^,v,k,(f)g,eQ,(i)Q,eK,(i)K} ■ (A2) 

These variables are related to the variables defined in Eq. ( p.2|) by 

^^^^ a^tan^^f^°^^V (A5) 



Q'x\ ^ ( cos (t)q 
qy) ^\sm(t)q 

kj [or 



(A6) 
(A7) 



There are only eight integrals because we find that upon changing to the above variables, 
the remaining integration can be done by symmetry. Here, T],(f)Q,(j)K, and 0g lie between 
and TT, while 9q and 9k lie between and tt/2. The variable ^ is used to set the upper limit 
for kz and qz- The scale factors aq and ax are taken to be 

ag = uj/-f ,aK = uj/-f . (A8) 

Monte Carlo methods reduce the Eq. { \i.'2\) to a summation 

= J2 fi^ii^^i2, Xis, g)A^x , (A9) 

where A"x is the volume element of the subregion. For a finite volume V, we divide the 
integral region to equal volume elements and obtain 

•^('?) " ^ ^ fixii,Xi2, Xis, q) (AlO) 

where N is the number of subvolume elements. We can also monitor the fluctuation of the 
results, and according to the central limit theorem for large values of N, we can write the 
associated error 

S' = {<f>-<f >')/N , (All) 

where 

< f >= J^Yl ^*2, Xi8, q) 

< f >= J^Yl fi^iU^i2, ■■; Xi8, q) ■ (A12) 



13 



REFERENCES 



[1] C. Bottcher and M. R. Strayer Nucl. Instrum. Meth. B 31, 122 (1988). 
[2] C. Bottcher and M. R. Strayer, Phys. Rev. D 39, 1330 (1989). 

[3] C. Bottcher and M. R. Strayer, J. Phys. G: Nucl. Part. Phys. 16, 975 (1990). 
[4] Quark Matter 1990, Nucl. Phys. A 525, edited by J. P. Blaizot, et al.. 
[5] K. Kajantie, M. Kataja, L. McLerran, and P.V. Ruuskanen, Phys. Rev. D 34, 811 
(1986). 

[6] K. Kajantie, J. Kapusta, L. McLerran, and A. Mekjian, Phys. Rev. D 34, 2746 (1986). 

[7] H. Gould, in Atomic Theory Workshop on Relativistic and QED Effects in Heavy Atoms, 
Gaithcrsburg, Maryland, 1985, edited by Hugh P. Kelly and Yong-Ki Kim (AIP Conf. 
Proc. 136), (AIP, New York, 1985). 

[8] C. Bottcher and M. R. Strayer, in Physics of Strong Fields, Proceedings of the Interna- 
tional Advanced Course, Maratea, Italy, 1986, edited by W. Greiner, (NATO ASI Ser. 
B., Vol. 153) (Plenum, New York, 1987), page 629. 

[9] E. Teller, in Proceedings of the Ninth International Conference on the Application of 
Accelerators in Research and Industry, edited by J. L. Duggan (North-Holland, New 
York, 1986). 

[10] C. Bottcher and M. R. Strayer, in Proceedings of the Second Workshop on Experiments 
and Detectors for RHIC, Berkeley California, May 25-29, LBL-24604 (1987). 

[11] M. J. Rhoades- Brown, T. Ludlam, J.-S. Wu, C. Bottcher, and M. R. Strayer, 
Fourth Workshop on Experiments and Detectors for a Relativistic Heavy-ion Collider, 
Brookhaven National Laboratory Report BNL 52262, 325 (July 1990). 

[12] J.-S. Wu, C. Bottcher, and M. R. Strayer, Phys. Lett. B252, 37 (1990). 

[13] J. A. M. Vermaseren, Nucl. Phys. B 229, 347 (1983). 

[14] C. R. Vane, S. Datz, P. F. Dittner, H. F. Krause, C. Bottcher, M. R. Strayer, R. Schuch, 

H. Gao, and R. Button, Phys. Rev. Lett. 69, 1911 (1992). 
[15] R. Baur et al, Nucl. Phys. A 566, 87c (1994). 

[16] A. Belkacem, H. Gould, B. Feinberg, R. Bossingham, and W. E. Meyerhof, Phys. Rev. 
Lett. 71, 1514 (1993). 

[17] Experiment WA99, approved by CERN SPS Program Committee March, 1993. Measure- 
ments of Pair Production and Electron Capture from the Continuum in Heavy Particle 
Collisions, Oak Ridge National Laboratory, Manne Siegbahn Institute of Physics, Univ. 
of Aarhus, and Univ. of Lund. 

[18] H. Gould, private communication. 

[19] P. B. Eby, Phys. Rev. A 43, 2258 (1991). 

[20] K. Hencken, D. Trautmann, and G. Baur, Phys. Rev. A 49, 1584 (1994). 

[21] M. Fatyga, M. J. Rhoades- Brown, and M. Tannenbaum, Can RHIC be used to test 

QED?, Brookhaven National Laboratory Workshop Proceedings, Upton, New York, 

BNL-52247, April 1990. 
[22] G. Baur, Phys. Rev. D 41, 3535 (1990). 

[23] M. Fatyga and J. W. Norbury, Fourth Workshop on Experiments and Detectors for a 
Relativistic Heavy-ion Collider, Brookhaven National Laboratory Report BNL 52262, 
345 (July 1990). 

[24] G. Baur, Phys. Rev. A 42, 5736 (1990). 

[25] M. J. Rhoades-Brown and J. Weneser Phys. Rev. A 44, 330 (1991). 



14 



[26] C. Best, W. Greiner, and G. Soff, Phys. Rev. A 46, 261 (1992). 

[27] K. Hencken, D. Trautmann, and G. Baur, Institut fiir theoretische Physik der Univer- 

sitiit Basel preprint, (1994). 
[28] C. F. von Weizsacker, Z. Phys. 88, 612 (1934). 

[29] E. J. Williams, and K. Dan. Widensk, Selsk. Mat. Fys. Medd.13, 4 (1935). 
[30] S. J. Brodsky, T. Kinoshita, and H.Terazawa, Phys.Rev. D 4,1532 (1971). 
[31] F. E. Low, Phys. Rev. 120, 582 (1960). 
[32] R. H. Dalitz and D. R. Yennie, Phys. Rev. 105, 1598 (1957). 

[33] D. L. Olsen, B. L. Barman, D. E. Greiner, H. H. Heckman, P. J. Lindstrom, G. D. 
Westfall, and H. J. Crawford, Phys. Rev. C 24, 1529 (1981). 

[34] R. Anholt and H. Gould, in Advances in Atomic and Molecular Physics, edited by B. 

Bederson (Academic, New York, 1987). 
[35] A. Goldberg, Nucl. Phys. A240, 636 (1984). 

[36] G. Bauer and C. A. Bertulani, Phys. Rev. C 35, 836 (1987); Phys. Rep. 161, 299 (1988). 

[37] J.D. Jackson, Classical Electrodynamics, (New York, Wiley, 1975). 

[38] V. B. Bcrcstctskii, E. M. Lifshitz, and L. P. Pitacvskii, Quantum Electrodynamics, 

second edition, translated by J. B. Sykes and J. S. Bell, (Pergamon Press, New York, 

1979), page 438. 
[39] C. A. Bertulani and G. Baur, Phys. Rep. 163, 299 (1988). 
[40] M. Vidovic, M.Greiner, C. Best, and G. Soff, Phys. Rev. C 47, 2308 (1993). 
[41] C. Bottcher, M. R. Strayer, C. J. Albert, and D. J. Ernst, Phys. Lett. 237B, 175 (1990); 

C. J. Albert, D. J. Ernst, C. Bottcher, and M. R. Strayer, in Progress in Nuclear Physics, 

ed. W-Y. Pauchy Hwang, S-C. Lee, C-E. Lee, and D. J. Ernst (North HoUand, New 

York, 1991). 

[42] A. S. Umar, M. R. Strayer, D. J. Ernst, and K. R. S. Devi, Phys. Rev. C 30, 1934 
(1984). 



15 



FIGURES 

FIG. 1. Direct (a) and crossed (b) Feynman diagrams for pair production in a heavy-ion colli- 
sion. 



FIG. 2. Schematic diagram depicting a relativistic heavy-ion collision. 

FIG. 3. The function ^(q) versus q for the production of electron-positron pairs. The points 
are the results of the Monte Carlo calculation and the smooth curve is our fit to these points. 



FIG. 4. The differential cross section da/dh versus impact parameter b for electron-positron 
pair production. Solid lines: exact numerical result; dot-dashed lines: equivalent-photon approxi- 
mation (Ref. p|). 



FIG. 5. The integrated cross section as a function of bmax for electron-positron pair production. 



Solid lines:exact numerical result;dot-dashed lines: equivalent-photon approximation (Ref. |25|) 



FIG. 6. Probability distribution for A^-pair production as a function of impact parameter h 
for Pb -|- Pb collision. Solid lines: exact numerical result; dot-dashed lines: equivalent-photon 
approximation (Ref. [^5|). Total A^-pair cross sections are given in units of barns. 



FIG. 7. Probability distribution for A^-pair production as a function of impact parameter 
h for S + Au collision. Solid lines: exact numerical result; dot-dashed lines: equivalent-photon 
approximation (Ref. [^5|). Total A^-pair cross sections are given in units of barns. 
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